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ABSTRACT 



The Method of Partial Images has been successfully 
applied to electrostatic problems involving conductors on 
a dielectric substrate. This same method is investigated 
for its adaptation to dynamic problems'. A Green's Function 
is derived and applied to the problem of a wire dipole 
antenna on a dielectric substrate. The input admittance 
of the antenna is computed by the Method of Moments. 
Experimentally measured values of input admittance are 
compared with the theoretical values and the error is 
discussed . 
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I . INTRODUCTION 



If a monopole antenna is mounted on a perfectly 
conducting ground plane, it will reflect totally above the 
plane and will produce an image below the plane. This 
image will have an image current in the same direction 
as the monopole in order to satisfy boundary conditions 
at the ground plane. In essence, because of its image, 
the monopole will act as a dipole in a free-space. Hence, 
if the ground plane is removed, the image may be replaced 
by another physically real monopole thereby creating a dipole. 
In other words, conditions above the reflective plane remain 
unchanged if the plane is removed and a real object replaces 
the image. 

A. THE USE OF IMAGE THEORY 

With the monopole, the thickness of the ground plane 
is assumed to negligible. This results in only one image. 
However, if the ground plane is a dielectric with finite 
thickness, total reflection no longer occurs. Instead, 
images are produced from partial reflection at the surface, 
reflection within the dielectric and transmission through 
the dielectric. These images are used in the "method of 
partial images" [1]. When the dielectric region is removed, 
these images may be considered real charges and are summed 
as an infinite series. However, a finite number of them 
give good results. 
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B. A DYNAMIC SOLUTION 



The method of partial images has been successfully 
applied to static problems. The purpose of this research 
is to investigate its application in solving dynamic 
problems. A search for a Green's Function is conducted. 
Once found, it is utilized in an impedance equation. 

The equation generates an impedance matrix so that the 
"method of moments" can be used to solve for the current 
distribution or the input impedance of an antenna mounted 
on a dielectric substrate. 
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II. NATURE OP THE PROBLEM 



Since the electrostatic problem has been solved, a 
method for solving the dynamic problem is sought using 
image theory. The dynamic situation is chosen to be a 
half-wave dipole antenna on a dielectric substrate of 
arbitrary thickness and relative permittivity (e^). A 
wire constitutes the antenna whose radius and length are 
variable . 

A. ON THE IMAGE COEFFICIENT METHOD 

Consider one small segment of the wire as an element 
of charge. Assume the charge to be a short distance from 
the dielectric. The charge element will have many lines 




of flux in every radial direction. However, in considering 
a single line of flux, some of it penetrates the dielectric, 
while some of it is reflected at the interface. The amount 
reflected is proportional to the reflection coefficient (K). 
Similarly, the part penetrating the dielectric is propor- 
tional to (1-K). 

Naturally, many reflections and transmissions occur 
due to this single flux line. In turn, it produces many 
image points. The diagram serves to illustrate: 
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x nr m 




K ' ~ \ - 

60+61 l * ir 



(Figure 1 . ) 



Each of the three regions, I, II, III, may. have different 
dielectric constants. However, in this case, regions I and 
III have £ r = l (free space). Region II has an arbitrary . 
In solving for the Green's Function, each region has a 
different image representation. 

(a) . Image representation for left region - I. 

(b) . Image representation for center region - II. 

(c) . Image representation for right region - III. 
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B. SOLVING FOR THE GREEN'S FUNCTION 



(<0 

U) 

(«) 



Specifically, in the problem at hand, the wire actually 
lies on the dielectric at the interface of regions II and 
III. Region II is the desired region with which to work 
because the accumulation of charge is assumed to be greater 
on the dielectric side of the wire than on the free space 
side. The image representation for this region appears as 



ir 

— Y^-^A 

-p 3 o-p)^ -p(i-p)^ o-p)<j pVp)<| p 4 (i-p)<^ 

where the reflection coefficient is now p c. 1 ~ £l in order 

l + € r 

to avoid confusion with the propagation constant (k). 

Adding terms to simplify: 

(a) - F°r t : o-p) 
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In terms of a summation, the Green's Function for the 
dielectric region is expressed t ^ ^ 
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Yet, with a single wire, variation occurs in only the 



X-direction . 
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The dielectric Green's Function then becomes 



* » i 



where 






C. PHASE CONSIDERATION 

In formulating the Green's Function for the static 
problem, the phase of each image point need not be consid- 
ered. However, in the dynamic problem, image point phase 
must be examined. With the charge element a distance "a" 
from the dielectric in Figure 1, each image point has a 
phase term associated with it. This means that ^d could 
not be expressed as a simple summation. 

On the other hand, phase terms are eliminated if the 
charge element is adjacent to the dielectric. Since the 
wire is actually adjoined to the dielectric, every image 
point appears as though it were immersed in the dielectric 
region. For instance, consider image point "q^" which 
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results from image path Cl. This Cl comes from flux line A, 
which is reflected into B and then C. The image path length 
(Cl) is equal to the addition of dielectric path lengths A 
and B through simple geometry. In other words, the' image 
path "Cl" appears to be in the dielectric medium due to A 
and B. Then image point "q^" also appears to originate in 
the dielectric for the same reason. This is true for every 
image point. Hence, all image points appear to be in the 
same medium with equal phase. The Green's Function, then, 
retains its simple form. 

D. THE USE OF 

The Green's Function <+d ) of this section was derived 
in a manner similar to Silvester's solution [1] using the 
method of partial images. Silvester, of course, solved a 
static problem. Yet, it seems logical to assume that his 
method should be applicable to a dynamic problem through 
the reasons stated. 

The method of .moments is used to find a numerical 
solution to the input impedance of the antenna. This method 
utilizes the Green's Function in solving for the impedance 
matrix. 
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III. THE IMPEDANCE MATRIX 



In solving for the current distribution or the feed- 
point impedance of an antenna, the matrix approach in the 
method of moments is a good approximation. The key is to 
properly’ load the impedance matrix. Once this is accom- 
plished, assumed values of voltage fill the voltage matrix 
and the problem is solved. 



The impedance matrix is obtained through the impedance 
equation. The equation is derived by Harrington [2] and 
the details are given in Appendix A. The Green’s Function 
(t d ) from the previous section is employed in the equation. 



A. CALCULATION OF THE GREEN’S FUNCTION 

In order to express the Green’s Function in matrix form, 
ij/ - ^ must be integrated along the wire. Holding the source 
point (n) fixed while the observation point (m) ranges over 
every segment generates the n— column of the matrix. Then 



[I] = [Z]" 1 [V] 




^(n,m) - y(n,m) + ^(n,rrv)J 
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S "t 

moving the source point to the (n+1) — segment and, again, 
allowing the observation point (m) to roam along the wire 
generates the next column, etc. Finally, the (m x n) 
matrix is filled, where "m" is chosen numerically equal 
to "n". 

In terms of equations : 
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There is a different Green's Function for each case. 



The derivations are detailed in Appendix B while the results 
are shown here. 

For case 1. (m=n) 
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For case II. (m^n) 
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B. CALCULATION OF THE IMPEDANCE MATRIX 

The impedance values are now found using the calculated 
Green's Functions in the impedance equation. 

Z ( m,n) = J Oil A$n (n ( m) + _J Fty (n m) 

** ^ (n | m) " m)J 
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Where A5 n ~ A$n’ since the length of each segment 

is chosen to be equal. 

Three different situations arise in calculating elements 
of the impedance matrix. The first is the "self-impedance" 
term. Whenever the observation point and source point 
coincide, a self-impedance term is generated. These terms 
occupy the diagonal of the matrix and are all numerically 
equal. 

The second situation is the "mutual-adjacent" term. 

These occur when the observation point is adjacent to the 
source point. Due to the fact that the distance between 
"m" and "n" is always the length of one segment ( AJU)> 
every mutual-adjacent term is equal. 

The third situation happens when the observation and 
source points are one or more segments apart and are called 
"mutual impedance" terms. They complete the matrix and 
vary in magnitude. 
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2. Mutual-Ad,- j acent Impedance 
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3. Mutual Impedance 
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C. PROGRAMMING 



The impedance equations of the previous section are 
programmed for computer usage. Of course, similarities 
exist throughout the equations making the programming 
easier. 

Once the proper parameters and constants are selected, 
the impedance matrix is filled and then inverted. The 
values of voltage in the voltage matrix are chosen to be 
zero everywhere except at the feedpoint, where the source 
is one volt. 

-1 

[Z] [V] = [I] 

The resulting matrix [I] will be the current distribution 
on the wire. The middle value of [I], the feedpoint value, 
is the feedpoint admittance since the source was one volt. 
The computer program is written to solve for various values 
of feedpoint, or input, admittances. 
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IV. DATA EVALUATION 



After programming the equations, meaningful data is 
sought for analysis. Several parameters are available: 
the radius, the length, the dielectric constant, the 
thickness, and the frequency. In order to evaluate the 
validity of the equations, a comparison is made against 
Harrington’s solutions.' 1 ' Hence, the dielectric constant 
in Region II is chosen to be e =1. The length is equal to 
half of a free space wavelength at f = 3 GHz. The length- 
to-diameter ratio, L/2a, is set at 7^.2, while the frequency 
is varied. 

The results are plotted and compare favorably to 
Harrington's graphs. Figures 3 and 4 show conductance and 
susceptance, respectively, versus the length-to-wavelength 
ratios. It was found that at least a forty-one by forty-one 
matrix was necessary to achieve accuracy. 

A. EXPERIMENTAL PROCEDURE 

Experimentally, a thin copper strip was attached to a 
one-eighth inch thick dielectric substrate with e =16. The 
dielectric was mounted perpendicular onto the ground plane 
so that the thin strip would be fed as a quarter-wave 
monopole vice a half-wave dipole. The antenna was driven 



"''Harrington, R.F., Field Computation by Moment Methods , 
p. 72, The MacMillan Company , 19687! 
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Susceptance (millimhos) Conductance (millimhos) 




(Figure 3.) 




Input Admittance of the Dipole 
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through 50 ohm miniature coax and input impedance was 
measured with a network analyzer. 

The experimental values of input conductance and sus- 
ceptance are plotted along with values calculated using the 
•impedance equations. Obviously, Figures 5 and 6 illustrate 
a significant difference between experimental and calculated 
results. There is an amplitude variation and a difference 
in resonant frequencies. 
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(Figure 5>) 
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V. CONCLUSIONS 



Prom the free-space curves which match Harrington's 
curves, it seems that the Green's Function and the 
impedance equations are valid. Of course, in this instance, 
e =1 so that the reflection coefficient (p) is zero and 
every term containing (p) is eliminated. The part remaining 
is simply a free-space Green's Pucntion. Hence, the solution 
is valid for a half-wave dipole antenna in free-space. 
Unfortunately, the same is not true when e >1. 



A. EXPERIMENTAL DIFFERENCES 

Aside from the normal attenuation’ losses of experimental 
measurements, it seems the only other inaccuracy would come 
from the adhesive used to attach the copper strip to the 
dielectric. This would influence the effective dielectric 
constant ( e gpp) and, perhaps, the amplitude of the admit- 
tance. Therefore, the emphasis is switched to the 
calculation procedures. 



B. PROBLEM IN THE GREEN'S FUNCTION 

The manner in which the Green's Function was chosen 
allows some error to exist. The accumulation of charge 
was assumed to be greater on the dielectric side of the 
wire rather than the free-space side. Because of this. 
Region II was selected in accordance with the method of 
partial images to solve for the necessary Green's Function. 
In the process, the propagation constant is dependent on 
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the relative permittivity of the dielectric. 



k= x 
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However, from the physical standpoint, as the limit of 
the dielectric thickness approaches zero, the Green's 



from section II-B will never reduce to a free-space function 
due to the dependence of the propagation constant on . 
This suggests the use of both a free-space propagation 
constant and a dielectric propagation constant. But there 
does not seem to be any apparent method in combining the two 
propagation constants with a Green's Function for Region II. 

Further investigation led to a Green's Function for 
Region III which incorporated both propagation constants. 

It reduces to the free-space function as zero dielectric 
thickness is approached. However, the solution leads to 
negative conductance. 

Finally, since the method of partial images can produce 
valid solutions for static and free-space dynamic problems, 
it would seem that with some additional changes this same 
method should solve all dynamic problems. Unfortunately, 
the final answer lies beyond the research efforts of this 
paper. The method of partial images, however, warrants the 
need for further investigation. 



Function should reduce to a free-space solution. 
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APPENDIX A 



DERIVATION OF THE IMPEDANCE EQUATION 

The impedance equation is developed in a manner similar 
to Harrington’s derivation. Formulating the problem: 
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The wire is broken into many segments subject to the 
following conditions: 

(1) Current flow is assumed to be along the wire axis. 

(2) To the axial component of E at the surface, apply 
boundary condition (A-5). 

With these assumptions, equations (A-l) through (A-4) 
become 
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where " " is the "length variable" along the axis of the 

wire . 
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Figure A-l illustrates the division of the wire into 

"n" segments. According to the notation, any n— segment 

_ + 
starts at n, is centered at n, and terminates at n. The 

distance between any n and n is denoted by , the 

length of the segment. Summing integrals over many small 

segments of AC n approximates an integral over the entire 

length. Equations (A -6) through ( A— 9 ) are then further 

approximated 



~ jU)Ajj(n') “ (m) - 1? (m) (A-10) 
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Similarly 
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substituting equation (A-13) into (A-12), and equation 
< A— 14 ) into an equation similar to (A-12) for $(m). 
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it is often convenient to express the integral parts of the 
above equations in the form of Green’s Functions. 
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Where Rm = the distance from the midpoint of the m-^ 
segment to the integration point on the n— - 
segment . 



Using this definition, equations (A-ll) , (A-15) and 
(A-l6) become 
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Now substituting equations (A-17), (A-l8) and (A-19) into 

equation (A-10) and noting that Aj(^)~ A(m) « A £ m 
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Two further approximations are made 

X (n-vO ~ Xcn)X(*\m) 
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With these two substituted, equation (A-20) is written as 
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or 

Eljj (rfl) A ft* * ^ ' Z(»v\n) X(n) 

where n ' 

2T(nfl ( rt) = jcoji A l 'A It« + _L [" ^(n,r?i) 

** ^n,^) “ ^ r (n, ri\) + i^(h,r^) J (A- 2A) 

Equation (A-2^) is the impedance equation which applies 
for the case of (m=n) self impedance, as well as for (m^n) 
mutual impedance. 

Since "m" roams along all "n" segments of the wire, 
there is a set of "n" linear equations. In matrix form: 
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or 

[V] = [Z] [I] 

If [V] and [Z] are known, the current distribution may be 
found, 

[I] = [Z]" 1 [V] 

as well as the feedpoint impedance. 



32 



APPENDIX B 



DERIVATION OF THE GREEN ’ S ' FUNCTION 



There is a distinct Green’s Function for each case: 
when the observation point coincides with the source point 
(m=n); and when the observation point and source do not 
coincide (m^n) . 

Case I. (m=n) 
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a. The first term must be integrated over one segment 




Using a Maclaurin series expansion on the first term: 
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since the first two terms of this Maclaurin expansion give 



reasonably good results. 
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b. Integration is now performed on the second term of 
the Green's Function. 
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Finally, the Green's Function for Case I, (m=n) is 
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APPENDIX D. COMPUTER PROGRAM 

COMPUTER APPLICATION ON THE METHOD 
OF IMAGES USING REGION II 



IMPLICIT R EA L*4 (K) 

IMPLICIT C0MPLEX*8 (C,Z) 

DIMENSION Z(41,41), ZMPI41), Y(41) 
R EAL*4 CAMPS ( 41 ) 

REAL V(41 )/20*0.0, 1.0,20*0.0/ 
COMPLEX CABS 

TFE THICKNESS (D) IN METERS 
D=G. 003175 

TFE DIELECTRIC CONSTANT (ER) 

ER= 16 . 0 

TFE RADIUS (A) IN METERS 
A =0 . COO 50 0 

20 READI5 , 30 » END=260 ) DATA 
3C FORMAT ( F5 . 3 ) 

NUM=4 1 

M I D= ( NU M+ 1 ) / 2 

THE LENGTH OF THE WIRE IN METERS 
ELNGTH=0. 13493750 

CLTA=ELNGTH/(FLOAT (NUM) ) 

F= ( (3E08 )*DATA)/ELNGTH 
P 1=3. 14159 
W = 2 . 0 * P I*F 
U = P I * 4 . 0*1 . E— 7 
E 0=1.0/ (36.0*PI ) *l.E-9 



EQUATION CONSTANTS 

RhC=( 1 . 0-ER ) / < 1 .0+ER) 

B = C . 0 

K = W*S QRT ( U *E0*ER ) 

A 1=W*U* DLT A**2 

A2=( 1.0-RH0)/(2.0*PI) 

A3= ( 1 . 0/DLTA ) *ALOG ( DL T A/ A ) 

A4=l,0/ ( 2 . 0* P I ) * ( 1 ,0-RH0)**2 
A5=2. 0/ ( W*E0*ER ) 

A6 = A2 / 2 .0 
A7=A4/2 . 0 
A 8 = A5 / 2 . 0 
A9 = K/2 . 0 
A10=DLTA**2 
A 1 1=4.0*0*42 
AK1 =K*D LT A 

CD1=CEXP(CMPLX(B,~AK1) J/DLTA 
AK2 = 2 .0 *AK 1 

CC2=CEXP(CMPLX( B,-AK2 ) )/ (2 . 0*DL T A ) 

C 1=CMPLX( B t A 1 ) 

C2=A3-CMPLX( B»A9) 

SUMMATION OF THE EXPONENTIAL WITH DELTA (DLTA) 

C3=CMPLX(B»B) 

R 2= 0 . 0 

DC 40 J=1 , 10 
J 1 = 2 . 0* J- 1 
J 2=J**2 
R1=RH0**J1 
R2=R2+R 1 
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X 1 = SQRT ( A LO+A L 1* J 2 ) 

X2=K*X1 

C4=CEXP(CMPLX(B,-X2) ) 

C5= ( R 1 *C4 ) /X 1 
C3=C3+C5 
4 C CONTINUE 
C 

c 

C SUMMATION OF THE EXPONENTIAL WITH 2*DELTA ( 2*DLT A ) 

C 

C6=CMPLX( B,B) 

R2 = 0. 0 

DO 50 1=1, 10 
I 1 = 2. 0*1-1 
I2=I**2 
R4=RH0**I 1 
R3=R3+R4 

X3=SQRT( 4. 0*A10+A1 1*12 ) 

X4=K*X3 

C8=CEXP(CM PLX(B» — X4 ) ) 

C9=(R4*C8) /X3 
C6=C6+C9 
50 CONTINUE 
C 
C 

C SUMMATION OF THE EXPONENTIAL WITH 2*I*D 

C 

C20=CMPLX( B,B) 

DC 60 1=1 v 10 
I 2 = 2. 0*1-1 
14= 1**2 
R10=RH0**I3 
X2C=2 .0*1*0 
X 2 1=K *X 20 

C21=CEX F ( CMPLX ( B,— X21 ) ) 

C 22 = ( R 1 0*C 2 1 ) /X20 
C 20 = C 20 +C2 2 
60 CONTINUE 
C 
C 

C CALCULATION OF THE IMPEDANCE MATRIX 

C 

DO 110 M=1,NUM 
DO 100 N=1 , NUM 
L= I AB S( M-N) 

IF(L.EQ.O) GO TO 90 
I F ( L . EQ . 1 ) GO TO 80 
C 

c 

c FOR MUTUAL IMPEDANCE 

C 

LP = L + 1 
' LM=L— 1 
AK3=K*L*DLTA 

CD3=CEXP( CMPLX ( B,-AK3) )/( L*DLTA ) 

AK4 = K*LP*DLT A 

C D4=C EX P( CMPLX (B,-AK4) ) / ( LP * DLT A ) 

AK5=K*LM*DLT A 

C D5=C EX P( CMPLX ( B,-AK5 ) )/( LM*DLT A ) 

C1C=CMPLX( B,B) 

C 13=CMPLX ( B, B) 

C 1 6=C MP LX ( B * B ) 

R 5= 0 . 0 

DC 70 J = 1 , 10 

J 11 = 2. 0*J- 1 

J 22=J**2 

R6=RH0**J1 1 

R 5=R 5+R 6 

A12=(L*DLTA)**2 

X5 = SQRT ( A12+A1 1* J22 ) 

X6=K*X5 

CU = CEXP(CMPLX(B,-X6) ) 
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C 12= ( R6*C1 1 ) / X5 
C 1 0=C 10+C12 
A13=( LP*DLTA)**2 
X7=SQRT ( A13+A11*J22 ) 

XE=K*X7 

C14=CEXP(CMPLX( B , -X8 ) ) 

C i 5= ( R 6*C 1 4 ) /X7 
C 13=C L3+C 15 
A14=( LM*DLTA )**2 
X 5 = SQ R T ( A14+A11*J22) 

X 10=K*X 9 

C17=CEX P( CMPLX( B,-X10 ) ) 

C 1 8= ( R6*C 17) /X9 
C 16= C 16 +C 1 8 
70 CONTINUE 

Z 2 1=01* ( A6*CD3-A7*C10) 

Z22=(CMPLX ( B, -A5 ) ) *( A6*CD3-A7*C 10) 
Z23=(CMPLX(B,A3) ) * ( A6*CD4-A7*C1 3 ) 
Z24=(CMPLX(B, A8) ) * ( A6*CD5-A7*C1 6 ) 
Z(M,N)= Z21 + Z22 + Z 23 + Z24 
GC TO 100 



FOR MUTUAL-ADJACENT IMPEDANCE 

80 Z31=C1* ( A6*CD1-A7*C3 ) 

Z32 = (CMPLX(8, -A5 ) ) * ( A6*CD1 -A7* C3 ) 
Z33 = ( CMPLX( B, A8) ) : $(A6*CD2-A7*C6 ) 
Z34 = ( C MPLX ( B , A8 ) ) * ( A2*C2-A7 *C20 ) 
Z(M,N) = Z31+Z32+Z33+Z34 
C-G TO 100 



FOR SELF IMPEDANCE 

90 Z 11=C1*(A2*C2-A7*C20) 

Z12=(CMPLX(B,-A5) ) * ( A 2*C 2- A 7*C 2 0 ) 
Z13=(CMPLX(B,A5 ) ) * ( A6 * CD 1- A7 *C3 ) 

Z ( M t N ) = Z11 + Z12+Z13 

100 CONTINUE 
110 CONTINUE 



THE MATRIX INVERSION 

CALL CM1N1 (NUM,Z, NUM, DETERM) 



MATRIX MULTIPLICATION OF Z-INVERSE AND V 

DC 130 1= L , NUM 
ZMP(I ) =(0.0, 0.0) 

DO 120 J = 1 , N UM 
Z M P ( I ) = ZMP ( I )+Z( I, J)*V(J) 

120 CONTINUE 
130 CONTINUE 

WRITE (6 ,140 ) 

14C FORMAT ( • 1* ) 

WRITE (6, 15 0) 

150 FOP MAT ( 1 5 ( / ) ) 

WRITEt 6,160) D , E R , A • 

160 FORMAT ( 10X, 'THE THICKNESS D • , 1 5 X , F 12 . 8 , / / , 1 OX , • THE 
1 ' DIELECTRIC CONSTANT E R • , 4X , F 1 2 . 8 , // , 10 X , « T HE RADIUS' 
2' A' , 18X,F 12. 8,//) 

W R I T E ( 6 , 170) W 

170 FCRMAT(//,10X,' THE RADIAN FREQUENCY W • , 10 X , E 1 3 . 7 , / / ) 

W R I T E ( 6 , 180) 

130 FORMAT ( 10X, • THE COMPLEX CURRENT AT THE FEEDPOINT IN • 
3' AMPS' , /) 

WRITE(6, 190) ZMP(MID) 

190 FORMAT (42X,EU .5 ,3X, El 1.5, 3X ) 

DO 200 1 = 1 , NUM 
CAMPS ( I ) = CABS ( Z M P ( 1 ) ) 
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200 CONTINUE 

WRITE (6, 210) 

210 FORMAT (//, 10X, 'THE MAGNITUDE OF THE CURRENT',/) 

WRITE (6,220) CAMPS(MID) 

220 FORMAT ( 42X » E 11. 5 ) 

WR ITE (6,230) DATA 

230 FORMAT! //,10X, 'FOR THE LENGTH/ L AMDA RATIO *,5X,F6.4,/} 
W R ITE ( 6 , 240 ) F 

240 FORMAT!//, 10X, 'THE FREQUENCY F • , 16X, IP E 1 C . 4 , / / ) 



CARD PRINT OUT - DATA FOR DRAW SUBROUTINE 
WRITE! 7 ,250 ZMP(MID) 

250 FORMAT! Ell. 5) 

GO TO 20 
260 CONTINUE 
STOP 
END 



SUBROUTINE CMIN1 ( N , A , ND I M , DETE RM ) 



CM IN 1 IS A SUBROUTINE WHICH WILL ACCEPT A SINGLE 
PRECISION COMPLEX MATRIX AND RETURNS THE INVERSE 
OF THE MATRIX IN ITS PLACE. 

N - THE ORDER OF THE MATRIX TO BE INVERTED 

A - COMPLEX SINGLE PRECISION INPUT MATRIX 

(DESTROYED). THE INVERSE OF A IS 
RETURNED IN ITS PLACE. 

ND I M - THE SIZE TO WHICH A IS DIMENSIONED 



COMPLEX A ( ND I M , N D I M ) ,PIVOT(100) , AM AX, T, SWAP, DETERM, U 
I NT EGER *4 IPIVOT! 100) , INDEX! 100 ,2) 

REAL TEMP,ALPHA(100) 

INITIALIZATION 

DETERM=CMPLX( 1.0, 0.0) 

DC 20 J = 1 , N 
ALPHA ( J ) = 0 . ODO 
DO 10 1 = 1, N 

10 ALPHA! J ) = ALPHA(J )+A( J, I ) *CONJ G ( A ( J , I ) } 

ALPHA! J) =SGRT( ALPHA! J ) ) 

20 IPIVOT (J)=0 
DC 600 1=1, N 

SEARCH FOR PIVOT ELEMENT 

AMAX=CMPLX!0. 0,0.0) 

DO 105 J = 1 , N 

IF ( IPIVOT (J )-l ) 60,105,60 
6 C DC 100 K = 1 , N 

IF ( IPI VOT(K )-l) 80, 100,740 
8 0 T EMP= AMAX*CON JG ( AMAX )-A(J,K)*CONJG(A(J,K)) 

I F ( TEMP ) 85 , 85, 100 

8 5 I R 0 W = J 
I CO LUM= K 
AMA X=A ( J , K ) 

100 CONTINUE 
105 CONTINUE 

' IPIVOT! ICOLUM) = IP I VOT ( ICOLUMJ + l 

INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL 

IF! IROW-ICOLUM) 140, 260, 140 

1 4 C DE TERM=DE TERM 
DO 200 L=1,N 
SWAP=A ( I ROW , L ) 
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A ( I ROW, L) = A( ICOLUM * L) 

200 A ( ICOLUM, L )=SWAP 
S k A P = A L PH A ( I ROW ) 

ALPHA ( I RO W ) =A LPH A ( ICOLUM) 

ALPHA ( ICOLUM )=SWAP 
260 INDEX (1,1) = I R OW 

INDEX! I , 2)=IC0LUM 
P IVOT (I ) = A( ICOLUM, ICOLUM ) 

U= PIVOT (I) 

TEMP=P I VOT ( I ) *CON JG ( PI VOT( I ) ) 

IF(TEMP) 330, 720, 330 
C 

C DIVIDE PIVOT ROW BY PIVOT ELEMENT 

C 

330 A ( ICOLUM, ICOLUM) =CMPLX(1. 0,0.0) 

DO 350 L=1,N 
U = P IVOT ( I ) 

350 A ( ICOLUM, L)=A( ICOLUM, L) /U 
C 

C REDUCE NON-PIVOT ROWS 

C 

380 DO 550 L 1= I , N 

I F ( L I- ICOLUM ) 400, 550, 400 
400 T=A( LI , ICOLUM) 

A (LI, ICOLUM) = CMP LX (C. 0,0.0) 

DC 450 L= 1 » N 
U=A( ICOLUM, L) 

450 A ( L 1 , L )=A(L1»L)-U*T 
550 CONTINUE 
60C CONTINUE 
C 

C INTERCHANGE COLUMNS 

C 

620 DO 710 1=1, N 
L = N + 1 - 1 

I F( INDEX! L , L) -INDEX(L,2) ) 630, 710, 630 
63C JRCW= I NDEX ( L , 1) 

JCCLUM= I NDEX ( L , 2 ) 

DO 705 K= 1 » N 
S W A P= A ( K, JROW ) 

A! K,JROW) =A(K, JCOLUM) 

A ( K, JCOLUM )=SWAP 
705 CONTINUE 
7 I C CONTINUE 
RETURN 

720 WRITE (6,730) 

73C FORMAT! 2 OH MATRIX IS SINGULAR) 

740 RETURN 
END 
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